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ABSTRACT 

Using high-resolution A-body simulations we study the stellar properties of dwarf spheroidal (dSph) 
galaxies resulting from the tidally induced morphological transformation of disky dwarfs on a cosmo- 
logically motivated eccentric orbit around the Milky Way. The dwarf galaxy models initially consist 
of an exponential stellar disk embedded in an extended spherical dark matter halo. Depending on 
the initial orientation of the disk with respect to the orbital plane, different final configurations are 
obtained. The least evolved dwarf is triaxial and retains a significant amount of rotation. The more 
evolved dwarfs are prolate spheroids with little rotation. We show that in this scenario the final 
density distribution of stars can be approximated by a simple modification of the Plummer law. The 
kinematics of the dwarfs is significantly different depending on the line of sight which has important 
implications for mapping the observed stellar velocity dispersions of dwarfs to subhalo circular veloc- 
ities. When the dwarfs are observed along the long axis, the measured velocity dispersion is higher 
and decreases faster with radius. In the case where rotation is significant, when viewed perpendicular 
to the long axis, the effect of minor axis rotation is detected, as expected for triaxial systems. We 
model the velocity dispersion profiles and rotation curves of the dwarfs under the assumption of con- 
stant mass-to-light ratio by solving the Jeans equations for spherical and axisymmetric systems and 
adjusting different sets of free parameters, including the total mass. Wc find that the mass is typically 
overestimated when the dwarf is seen along the long axis and underestimated when the observation is 
along the short or intermediate axis. For the studied cases the effect of non-sphericity cannot however 
bias the inferred mass by more than 60 percent in either direction, even for the most strongly stripped 
dwarf which is close to disruption. 

Subject headings: galaxies: dwarf - galaxies: fundamental parameters - galaxies: kinematics and 
dynamics - galaxies: Local Group - galaxies: structure - cosmology: dark matter 



1. INTRODUCTION 

Dwarf spheroidal (dSph) galaxies of the Local Group 
(for a review see Mateo 1998) can provide crucial tests of 
the currently favored cold dark matter (CDM) paradigm 
of cosmological structure formation (e.g., White & Rees 
1978; Blumenthal et al. 1984). Several scenarios have 
been proposed so far for the origin of these systems, in- 
cluding gravitational interactions in the form of galaxy 
harassment (e.g., Moore et al. 1996) tidal stirring (Mayer 
et al. 2001; Kravtsov et al. 2004; Kazantzidis et al. 
2004b) or resonant stripping (D'Onghia et al. 2009) as 
well as hydrodynamical processes (e.g., Ricotti & Gnedin 
2005; Mayer et al. 2007; Tassis et al. 2008). The tidal 
stirring scenario proposed by Mayer et al. (2001, see also 
Gnedin et al. 1999; Peharrubia et al. 2008) can explain 
the morphological transformation of disky dwarf galax- 
ies resembling dwarf irregulars into pressure-supported 
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stellar systems with the structural properties of dSphs. 
When combined with ram pressure stripping and the ef- 
fect of the UV background (Bullock et al. 2000; Susa 
& Umemura 2004) it can also account for the low gas 
content and the very high mass-to-light ratios (M/L) of 
some of them (Mayer et al. 2006; 2007). 

In order to explore this scenario in more detail Kli- 
mentowski et al. (2007, 2009a) have recently studied the 
evolution of two-component dwarf galaxies orbiting the 
Milky Way on a cosmologically motivated eccentric orbit 
using high resolution A-body simulations. The dwarfs 
were initially composed of a stellar disk embedded in 
a more extended dark matter halo. Three simulations 
were performed which differed by the initial inclination 
of the disk with respect to the orbital plane. Although 
the evolution proceeded a little differently in each case, 
strong mass loss and morphological transformation of the 
disk always took place. After the first or second (out 
of five) pericenter passages the disk evolves into a bar 
which in the subsequent evolution becomes shorter. The 
morphological transformation of the stellar component is 
accompanied by the loss of angular momentum so that 
little rotation remains at the end. The final products are 
very similar in shape to the dSph galaxies of the Local 
Group. In this paper we study in detail the structural 
and kinematic properties of stars in these dwarfs at the 
end of their tidal evolution. 

DSph galaxies are also interesting because of their sig- 
nificant dark matter content. It turns out that the tidal 
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TABLE 1 

Properties of the simulated dwarfs. 



simulation 


i[deg] 


»*max [kpc] 


M tot (r max )[M ] 


N stars ( ^ max ) 


ArjM^max) 


M stars (r max )[M ] 


M DM (r ma x)[M ] 


A 


90 


2.00 


3.73 x 10 7 


73423 


25394 


1.10 x 10 7 


2.63 x 10 7 


B 


45 


1.70 


2.05 x 10 7 


37739 


14378 


5.64 x 10 s 


1.49 x 10 7 


C 





1.44 


4.67 x 10 6 


8379 


3298 


1.25 x 10 6 


3.42 x 10 6 



stirring scenario is able to reproduce both the mildly dark 
matter dominated dwarfs like Fornax or Leo I (Klimen- 
towski et al. 2007, 2009a) and the very dark ones like 
Draco (Mayer et al. 2007). In spite of significant mass 
loss due to tidal stripping the dwarfs are able to retain 
enough dark matter to remain mostly dark, although the 
extended dark halo is lost very fast and the mass approx- 
imately traces light in the final stages. This means that 
according to this scenario the mass modeling of dSph 
galaxies can be reliably done adopting the scale lengths 
of the distribution of light as is indeed a common ap- 
proach. 

The presence of tidal stripping however also means that 
in addition to the contamination from Milky Way stars, 
which can be dealt with using photometric methods (Ma- 
jewski et al. 2000) or metallicities (Walker et al. 2009a), 
the kinematic samples used for dynamical modeling can 
be contaminated by tidal debris. Klimentowski et al. 
(2009b) have shown that such contamination is indeed 
very probable because of the typically radial orientation 
of tidal tails for dwarfs on eccentric orbits. This contami- 
nation may artificially inflate the velocity dispersion and 
bias the inferred parameters of the density profile or the 
anisotropy of stellar orbits. It can however be reliably 
removed using procedures for interloper rejection (Kli- 
mentowski et al. 2007; Lokas et al. 2008). 

Another source of systematic error in the dynamical 
modeling may be due to departures from sphericity. As 
mentioned above, the products of tidal stirring are not 
spherical but rather form prolate spheroids. In this 
work we study the effect of non-sphericity on mass and 
anisotropy estimates, i.e. we show what biases are ex- 
pected if spherical models are applied to galaxies whose 
shapes depart from spherical. We also discuss how the 
estimates can improve if axisymmetric models are used 
in those cases where rotation is detected. 

The paper is organized as follows. In section 2 we pro- 
vide a brief description of the simulations used in this 
work. In section 3 we study the properties of the stel- 
lar component of three dSph galaxy models obtained in 
three simulations with different initial inclination of the 
disk. Section 4 discusses the line-of-sight density and ve- 
locity distributions, as they would be seen by a distant 
observer. In section 5 we briefly compare the proper- 
ties of our simulated dwarfs to observations. Section 6 
is devoted to modeling the mock kinematic data using 
spherical and axisymmetric models and studying the bi- 
ases inherent in such modeling for different lines of sight. 
The discussion follows in section 7. 

2. THE SIMULATIONS 

The simulations used for the present study follow the 
evolution of dwarf disk-like systems orbiting within the 
static host potential of a Milky Way-sized galaxy and are 
described in detail in Klimentowski et al. (2007, 2009a). 
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Fig. 1. — Density distributions of stars in the simulated dwarfs. 
The first, second and third row corresponds respectively to dwarf 
A, B and C of Table [T] The left column panels show the density 
distribution of stars measured in radial bins (filled circles) and 
the solid lines show the Plummer law JTJ with parameters from 
Tablc[2]fitted to the data. The second column presents the densities 
measured in narrow cuboids along the x (circles), y (triangles) and 
z (squares) axis of the dwarfs. The solid, dashed and dotted lines 
show respectively the best-fitting profiles obtained with density 
following from the potential with parameters listed in Table [3] 



Here we provide a short summary of the most important 
parameters. Live dwarf galaxy models are constructed 
using the technique by Hcrnquist (1993) and consist of 
an exponential stellar disk embedded in a spherical and 
isotropic Navarro et al. (1996, hereafter NFW) dark mat- 
ter halo. The structural properties of the dark halo and 
disk are related through disk galaxy formation models 
in the currently favored concordance CDM cosmological 
model (Mo et al. 1998). 

The dwarf progenitor has the total mass M — 4.3 x 10 9 
Mq. The mass and radial scale length of the disk are 
1.5 x 10 8 M Q and i?d = 1.3 kpc, respectively, and its 
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TABLE 2 

Fitted parameters of the spherical models for the 
distribution of stars. 



TABLE 3 

Fitted parameters of the non-spherical models for the 
distribution of stars. 



simulation 


A4[M ] 


a [kpc] 


simulation 


M B [M Q ] 


a [kpc] 


b [kpc] 


c [kpc] 


d [kpc] 


A 


1.33 x 10 7 


0.73 


A 


1.62 x 10 7 


0.83 


0.93 


0.81 


1.41 


B 


7.01 x 10 6 


0.66 


B 


8.74 x 10 6 


0.76 


0.74 


1.16 


1.16 


C 


2.00 x 10 6 


0.86 


C 


2.70 x 10 6 


1.02 


0.77 


1.22 


1.22 
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Fig. 2. — Left panel: the final mass profiles of the three simulated 
dwarfs A, B and C (solid, dashed and dotted line respectively). 
The vertical dashed line indicates the values of mass within 300 
pc, M3oo- Right panel: the final mass-to-light density ratios for 
the three dwarfs. 

vertical structure is modeled by isothermal sheets (with 
the vertical scale height z& — 0.13 kpc). The halo virial 
mass and concentration are M V j r = 3.7 x 10 9 M Q and 
c = 15, respectively. The halo is exponentially trun- 
cated outside the virial radius r v i r — 40.2 kpc to keep 
the total mass finite (Kazantzidis et al. 2004a) and adia- 
batically contracted in response to the growth of the disk 
(Blumcnthal et al. 1986). The host galaxy is modeled 
by a static gravitational potential based on the dynam- 
ical mass model Al for the Milky Way (Klypin et al. 
2002). It consists of a NFW halo with the virial mass of 
.Mvir = 10 12 M Q and concentration c = 12, a stellar disk 
with mass M D = 4 x 10 10 M , the scale length R d = 3.5 
kpc and the scale height z d = 0.35 kpc, and a bulge of 
mass Mb = 0.008Af vir and scale-length a h = 0.2R d . The 
dwarf galaxy evolves on an eccentric orbit with apoc- 
enter r a = 120 kpc and pericenter to apocenter ratio 
of r p /r a rs 0.2, close to the median ratio of pericentric 
to apocentric radii found in high-resolution cosmological 
./V-body simulations (Ghigna et al. 1998; Diemand et al. 
2007). The evolution is followed for 10 Gyr correspond- 
ing to approximately five orbital times. We used three 
simulations with the disk initially inclined by 0°, 45°, 
and 90° with respect to the orbital plane. 

Lastly, the simulations were performed using PKD- 
GRAV, a multistepping, parallel, tree TV-body code 
(Stadel 2001). We sampled the live dwarf galaxy with 
4 x 10 6 dark matter particles and 10 6 stellar particles, 
and employed a gravitational softening length of 100 and 
50 pc, respectively. 

3. THE PRODUCTS OF TIDAL STIRRING 

Table [T] lists the basic properties of the simulated 
dwarfs in the final stage. For a detailed description of 
the intermediate evolution we refer the reader to Kli- 
mentowski et al. (2009a). In order to avoid any contam- 
ination from the tidally disrupted stars we select only the 



stellar and dark matter particles within radius r < r max 
kpc from the center of the dwarf. The adopted values 
of r mal , listed in the Table, were determined by inspec- 
tion of the stellar density profiles, i.e. they are the radii 
at which the profiles start to flatten signifying the tran- 
sition to the tidal tails. The radii are of the order of 
the dwarf's tidal radius or slightly smaller. The Table 
also gives the initial inclinations of the stellar disk plane 
with respect to the orbital plane, the total masses of the 
dwarfs within r max , the numbers of stars and dark parti- 
cles contained inside these radii and the stellar and dark 
masses separately. In section 5 we discuss the observa- 
tional properties of our dwarfs. 

The left column of Figure [T] shows with the filled circles 
the radially averaged density profile of the stars in the 
three simulated dwarfs measured in the range between 
r = 0.2 kpc (which is larger than our resolution limit) 
and r = r max in logarithmic bins. The data are well 
fitted by the Plummer distribution 



p{r) 



3a 2 M s 



1 



4tt (r 2 + a 2 ) 5 / 2 ' 



(1) 



where M s is the total mass of the stellar component and 
a is the characteristic scale-length. This density follows 
via the Poisson equation from the potential 

•M-^fV P) 

The fits were done by minimizing the sum Ei(pdata — 
AnodcOVPmodei- ^ e best-fitting profiles are plotted as 
solid lines in the left-column panels of Figure Q] and the 
best-fitting parameters are listed in Tabled The corre- 
sponding half-light radii are ri/2 = 1,0.9 and 1.15 kpc 
respectively for model A, B and C. 

In order to describe the departures from sphericity in 
the structure of the dwarfs using the stellar particles we 
determine the principal axes of the inertia tensor and 
define our coordinate system so that the x axis lies along 
the longest axis of the stellar distribution and the z axis 
along the shortest. We measure the density along the 
X, y and z axis in cuboids of logarithmic bins along the 
given axis and of ±0.2 kpc width along the remaining 
dimensions, so for example for the measurements along 
the z axis we adopt |a;| < 0.2 kpc and \y\ < 0.2 kpc. The 
measurements were made along the positive and negative 
side of a given axis and then averaged. 

The density profiles measured along the three axes are 
shown in the right column of Figure [T] with circles corre- 
sponding to the measurements along the x axis, triangles 
along the y axis and squares along z. We see that the 
shape of model A is clearly triaxial: the densities are 
significantly different along the three axes: the distribu- 
tion is rather flat along the x axis signifying the presence 
of a remnant bar and very steep (almost a power law) 
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Fig. 3. — Mean velocities along the three axes v x (dotted line), v y (solid line) and v z (dashed line) as a function of x, y and z (left to 
right column). The three rows show results for dwarf A, B and C respectively. 



along the z axis. As we demonstrate below, the z axis 
almost coincides with the rotation axis of the dwarf. This 
means that in addition to the presence of the remnant 
bar, dwarf A is also flattened due to rotation. 

The shapes of the models B and C are spheroidal, i.e. 
the density distribution is rather flat along the x axis 
and steeper along y and z but the latter two profiles are 
indistinguishable (therefore our choice of y and z is ran- 
dom with respect to x). Note also that the distribution 
of stars in model C is noisier, the dwarf is irregular and 
close to disruption, as indicated by its much smaller mass 
(see Table [I]). 

The density distribution along the three axes can be 
well approximated by the density following by the Pois- 
son equation from the modified Plummer potential 

<Z>(x,y,z) = (3) 

GMs 

[(a; 2 + y 2 + z 2 + a 2 ) 2 - b 2 x 2 + c 2 y 2 + d 2 z 2 } 1 / 4 

where again M s is the total mass while a, 6, c and d are 
constants with the dimension of length. The first term in 
the denominator is the same as in the original Plummer 



potential, (r 



-). The second, negative term, —b 



introduces a bar-like distribution along the x axis and 
the third and fourth positive terms c 2 y 2 and d 2 z 2 , make 
the profile steeper along the y and z axis respectively. 
The form of the potential is a generalization of the 
modified Plummer distribution proposed by Lynden-Bell 



(1962). 

We fitted the profiles measured from the simulation 
with the density following from ([3]) via the Poisson equa- 
tion by the same minimization scheme, except that now 
the data points for model A were weighted by the density 
measured along the y axis which have intermediate val- 
ues (otherwise the measurements along the z axis would 
be fitted best). For models B and C the data points 
were weighted by the average value of density along x 
and y. The solid, dashed and dotted lines in the second 
column of Figure[T]show respectively the best-fitting pro- 
files along the x, y and z axis respectively, obtained with 
parameters listed in Tabled Note that for models B and 
C we kept c = d. 

The important dynamical properties of the three sim- 
ulated dwarfs are illustrated further in Figure [2] This 
figure shows the final mass, M(< r), and mass-to-light 
density, p/l, profiles plotted up to each corresponding 
'"max- It can immediately be seen that the mass pro- 
files are very different in spite of the fact that the three 
dwarfs had identical initial density distributions and or- 
bits and differed only by the initial inclination of the disk 
with respect to the orbital plane. Interestingly, while our 
dwarfs started with identical masses within 300 pc, M300, 
they become substantially different among the three fi- 
nal models being almost an order of magnitude larger in 
dwarf A compared to that of dwarf C. 

Figure [5] also shows that the final mass-to-light den- 
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Fig. 4. — Dispersions of the velocities along the three axes <r x (dotted line), cr y (solid line) and cr z (dashed line) as a function of x, y and 
z (left to right column). The three rows show results for dwarf A, B and C respectively. 



sity ratios are not significantly different among the three 
dwarfs and, in addition, they are all almost constant with 
radius. The former is related to the fact that the initial 
mass-to-light distributions were identical and indicates 
that the dark matter and stars are lost in proportional 
amounts in all cases. The latter is a consequence of 
two facts. First, the initial mass-to-light distributions 
were constant within ~ 2 kpc or so (Klimentowski et al. 
2007). Second, after the outer part of the dark matter 
halo, which was initially much more extended than the 
stars, is stripped at the first pericenter passage, the stars 
and dark matter are stripped from all radii. Note that 
a similar result is obtained in simulations including gas 
dynamics (Mayer et al. 2007). The final mass-to- light 
profile could however be different if other initial condi- 
tions were assumed (e.g. Peharrubia et al. 2008). 

The velocity distribution of the stellar components of 
the dwarfs is illustrated in Figures [3] and [4] In Figure [3] 
we plot the mean velocities along the three axes, v^, 
and as a function of x, y and z. As in the case of 
density, the measurements were made in narrow cuboids 
placed along a given axis, e.g. for the measurements as 
a function of z we selected stars with \x\ < 0.2 kpc and 
\y\ < 0.2 kpc. For model A significant rotation can be 
seen, but mostly in the xy plane, i.e. much less rotation 
is seen along z. For model B rotation is very small and 



for model C negligible. This can be understood as due to 
more effective removal of angular momentum for models 
where the stars initially had angular momentum more 
aligned with the orbital angular momentum of the dwarf 
(e.g., Read et al. 2006). 

Figure0]shows the velocity dispersion profiles along the 
three axes a x , a v and a z as a function of x, y and z. Here 
and below the dispersions are measured with respect to 
the mean velocity in a given bin and calculated using 
the unbiased estimator s = EIL=i( u « — v) 2 /(n — I)] 1 / 2 
where n is the number of stars in a bin. We can see that 
the velocity dispersions along x are always largest, while 
those along y and z are smaller and comparable to each 
other. 

4. THE LINE-OF-SIGHT DISTRIBUTIONS 

In this section we discuss the line-of-sight distributions 
of density and velocity of the stars which can be directly 
compared to observations. The triaxiality of dwarf A is 
further illustrated in Figure [5] where we plot the surface 
density distributions of the stars seen along the three 
axes x, y and z from the top to the bottom row. The 
left column panels show the isodensity contours of the 
surface density measured from the simulation and the 
right column panels the corresponding values from the 
best-fitting model following from the modified Plummer 
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Fig. 5. — Surface mass distribution of the stellar component 
of dwarf A observed along the x, y and z axis (from the upper 
to the lower row). Left panels show the contours of equal surface 
density of stars measured from the simulation. Right panels plot 
the corresponding distributions of the modified Plummer model (J3j 
with parameters from Table [3] adjusted to the density profiles in 
Figure[T] The surface densities £ were expressed in Mq kpc -2 and 
the contours are spaced by AlogS = 0.2. The innermost contour 
level is logS = 7 (upper row), 6.8 (middle row) and 6.6 (lower 
row). 



potential ^ with parameters from Table [31 We see that 
the model reproduces quite well the basic features of the 
simulated dwarf. When viewed along the z axis (lower 
row) the remnant bar is well visible but there is otherwise 
little flattening. When viewed along the y axis (middle 
row), both the bar and the flattening along z are well 
visible. When viewed along the x axis (upper row) the 
bar is not seen, the contours are almost circular except 
for some flattening along the z axis. 

Figures [5] and [7] present analogous results for models 
B and C, except that now we show only the view along 
the longest axis (upper row) and along one of the other 
two axes (lower row) since the two look very similar. 
As expected for spheroidal shapes, the surface density 
contours for the view along the x axis are circular. In 
the perpendicular direction an elongated shape is seen. 
The two models look quite similar, although model B is 
more massive (see the different contour levels), except 
for the fact that in the case of model C the contours are 
much more irregular, which is a sign of a stage close to 
disruption. 

Figure [5] illustrates the line-of-sight kinematics of the 
stars in the simulated dwarfs. The measurements are 
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Fig. 6. — Surface mass distribution of the stellar component of 
dwarf B observed along the x axis (upper row) and perpendicular 
to it (lower row). Left panels show the contours of equal surface 
density of stars measured from the simulation. Right panels plot 
the corresponding distributions of the modified Plummer model J3} 
with parameters from Table [3] adjusted to the density profiles in 
Figure[T] The surface densities £ were expressed in Mq kpc~ 2 and 
the contours are spaced by AlogS = 0.2. The innermost contour 
level is logS = 6.8 (upper row) and 6.4 (lower row). 
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Fig. 7. — The same as Figure \E\ but for dwarf C. The contour 
spacing is A log S = 0.1 and the innermost contour level is log £ = 
6 (upper row) and 5.8 (lower row). 

done in the same configuration as for the surface density 
in Figures EE Note that no clipping of the data is nec- 
essary since we still consider only stars within r max . The 
columns from the left to the right correspond to observa- 
tions along the x, y and z axis. The first row shows the 
contours of equal line-of-sight velocity and the second 
row the contours of equal line-of-sight velocity disper- 
sion for model A. Next rows present analogous measure- 
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Fig. 8. — Line-of-sight kinematics of the stellar component observed along the x, y and z axis (from the left to the right column). The 
first and second row present the mean velocity and velocity dispersion maps respectively, for the dwarf in simulation A. The subsequent 
rows show the results for simulations B and C as marked on the right hand side of the Figure. Contour levels in km s -1 are indicated in 
each panel. Contours of negative mean velocity are shown with dashed lines. The measurements were made by binning the data in cells of 
size 0.28, 0.3 and 0.33 kpc respectively for simulation A, B and C respectively. 



ments for models B and C. Contour levels of the quan- 
tities expressed in km s _1 are marked in each panel and 
the dashed contour lines indicate negative values. 

Again, for model A rotation is well visible when the 
observation is along the x axis (left column) and along 
y axis (middle column) while it is much smaller for the 
line of sight along the z axis (right column). The velocity 
dispersions are similar for the observations along z and 
y while the dispersion measured along the longest axis is 
larger and decreases more steeply with distance from the 



center of the dwarf. A similar steeper decline of velocity 
dispersion along the longest axis is seen for model B, but 
the rotation level is low in this case. For model C the 
dispersion is rather flat even along the longest axis and 
rotation is very low. 

The velocity structure seems to be the richest in the 
case of model A, as expected for triaxial systems. Fig- 
ures [3] and [8] show that the orbital structure of our simu- 
lated dwarfs is similar to the one found for simple triax- 
ial models corresponding to Stackel potentials (de Zeeuw 
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TABLE 4 

Observational properties of the simulated dwarfs. 



simulation 


M v 


Hy [mag/arcsec 2 ] 


r i/2 [kpc] 


<t o [km s - 1 ] 


V/er 


A 


-11.8 


25.8-24.7 


0.53-0.88 


5-6.5 


0.16-0.62 


B 


-11.1 


26.0-25.2 


0.47-0.74 


4.5-5.5 


0.09-0.18 


C 


-9.7 


27.9-27.4 


0.73-0.98 


1.9-2.8 


0.11-0.25 



1985; Statler 1987). For such models four families of or- 
bits are identified: box orbits, short-axis tube orbits and 
inner and outer long-axis tube orbits. The only orbits 
with net angular momentum circulate cither around the 
short or the long axis in agreement with the rotation 
curves shown in Figure [3] The total angular momentum 
is thus generally misaligned from the short axis which re- 
sults in a velocity gradient along the apparent minor axis 
when the galaxy is observed in projection. This means 
that the zero velocity curve is misaligned with the minor 
axis, as we indeed see in the middle panel of the upper 
row of Figure [8] (observation along y) . The rotation is 
however much stronger around the short axis. For pro- 
late spheroids (models B and C) only the long-axis tubes 
are allowed but since the rotation levels are low in these 
cases the stars must travel in both directions so that no 
single (clockwise or anti-clockwise) direction dominates. 

5. COMPARISON WITH OBSERVATIONS 

To summarize the properties of the simulated dwarfs 
described in the two previous sections, in Table[4]we pro- 
vide the basic parameters that can be directly compared 
to observations. The visual magnitudes My were calcu- 
lated from the total stellar mass of the models assuming 
the stellar M/Ly = 3 M Q /L Q . For the remaining quan- 
tities, the central surface brightness fly, the half-light 
radius r 1 / 2 , central velocity dispersion uo and the ratio 
of the rotation velocity to the velocity dispersion V/ a we 
provide ranges of values corresponding to different lines 
of sight (along x, y and z). The values of V/a were calcu- 
lated by taking the maximum line-of-sight rotation and 
dispersion measured along the major or minor axis of the 
galaxy, as in Figure [5J 

Comparison with the observed values of these param- 
eters for the Local Group dSph galaxies (e.g. Mateo 
1998; Grebel et al. 2003; Tolstoy et al. 2009) shows that 
they are rather representative for this type of objects. 
Only the velocity dispersion values may appear a little 
too low. Note however, that we have considered here 
the final stages of the evolution, i.e. after 10 Gyr when 
the dwarfs have lost majority of their mass. Although 
the progenitors of dSph galaxies are expected to become 
satellites early, not all such galaxies have been tidally 
stripped for so long. Moving back 2 or 4 Gyr in our sim- 
ulations we would easily find earlier stages with velocity 
dispersions 50-100 percent higher but other properties 
very similar (see section 4 and fig. 12 in Klimentowski 
et al. 2009a). 

The parameters in Table [4] also obey the relationships 
like My-[iy or My-r\i 2 , i.e. in plots showing these re- 
lations (such as fig. 34 in Kormendy et al. 2009 or 
fig. 1 in Tolstoy et al. 2009) our dwarfs would occupy 
regions characteristic of dSph galaxies. Although with 
just three cases considered we cannot really explore cor- 



relations between the parameters, we can immediately 
see that there is a clear trend of surface brightness de- 
creasing with decreasing luminosity as observed in real 
dSph galaxies, in contrast to ellipticals (fig. 1 in Kor- 
mendy et al. 2009). This feature, among others, has been 
interpreted as pointing towards different formation sce- 
narios of spheroidal and elliptical galaxies by Kormendy 
(1985) and Kormendy et al. (2009). They suggested 
that elliptical galaxies may form mostly via mergers while 
sphcroidals are rather late-type systems that underwent 
transformation due to star-formation processes or envi- 
ronmental effects such as tidal stirring. This idea has 
recently gained support from the theoretical side; based 
on the analysis of the simulated Local Group Klimen- 
towski et al. (2009c) found that mergers of subhaloes 
are quite rare in such systems and occur early on so they 
cannot significantly contribute to the formation of a large 
fraction of dSph galaxies. The tidal stirring scenario thus 
seems to be the most effective gravitational mechanism 
by which such objects could form. 

Although our initial conditions constitute just one pos- 
sible model for the progenitors of dSphs, our results sug- 
gest that tidal interactions with the primary alone may 
introduce an appreciable scatter in M 300 (see Figure [2]), 
which could lead to some tension with claims of a com- 
mon mass scale for Galactic satellites within 300 pc (Stri- 
gari et al. 2008). In addition, while we cannot exclude 
the effects of baryonic physics on our findings, one could 
expect that a wider spectrum of initial conditions and 
orbits would likely further differentiate the inner masses 
of dSph galaxies. 

Recently, Walker et al. (2009c) suggested instead the 
existence of a scaling relation for dSph galaxies between 
the half-light radius and the mass contained within it of 
the form M(r 1/2 )/M Q = 5800(r 1/2 /pc) 14 . The masses 
M(ri/ 2 ) of our dwarfs (in units of 10 7 M Q ) are within 
the ranges 0.94-2, 0.52-1 and 0.22-0.33 (depending on 
the line of sight) respectively for dwarfs A, B and C. 
Combining with the half-light radii from Table |4] we find 
that our dwarfs fall below the line proposed by Walker et 
al., i.e. for these half-light radii they have smaller masses 
than the relation would predict. Note however that the 
masses used by Walker et al. may be overestimated due 
to contamination. With just three cases considered here 
(including one close to disruption) we cannot fully ad- 
dress the issue and we defer such investigations to future 
work. 

6. THE EFFECT OF NON-SPHERICITY ON MASS 
ESTIMATES 

6.1. Velocity moments and the Jeans equations 

In order to further study the observational conse- 
quences of non-sphericity of dSph galaxies we now at- 
tempt to model the kinematics of our simulated dwarfs. 
For each of the characteristic lines of sight discussed 
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Fig. 9. — Line-of-sight velocity dispersion profiles as a function of the projected radius R for dwarfs obtained in simulation A, B and C 
(from the upper to the lower row) for the observation along the x, y and z axis (from the left to the right column). The dotted lines show 
the best fitting isotropic (/3 = 0) solutions of the Jeans equation for the one-component model with scale length a adopted from the fit of 
the surface density of stars (Table [5] and Figure IT2l l and only the mass M adjusted. The solid lines are for models with a assumed and M 
and anisotropy /3 fitted. The dashed lines show the best-fitting isotropic models with M and a adjusted. In the upper left panel the dotted 
line coincides with the solid one. The best- fitting parameters are listed in Table [6] 



above, i.e. along the x, y and z axis, we measure the ro- 
tation velocity V and the velocity dispersion a along the 
line of sight. For this purpose we use all the available 
stars although in real observations the presently avail- 
able samples are at least an order of magnitude less nu- 
merous. We are however more interested here in study- 
ing systematic errors in the estimated parameters rather 
than estimating their realistic statistical errors. We as- 
sign to the measured values of the velocity moments the 
true sampling errors which are used to weigh the data 
points when fitting. Although these errors are unrealis- 
tically small in our case, for real samples they would be 
relevant, i.e. they would be larger by a constant factor 
and data points would be weighted by the same relative 
weights. 

In the cases where the rotation level is very low we pro- 
ceed in the standard way: we bin the stars in ten equal- 
number bins along the projected radius R (the distance 
from the center of the dwarf to the star on the surface 
perpendicular to the line of sight). The velocity disper- 
sion in each bin is calculated using the standard unbiased 
estimator of dispersion and assigned a sampling error of 
size a J \J1(n — 1) where n is the number of stars per bin. 



In the cases when rotation is detected (observation along 
the x and y axis for model A) we first transform the pro- 
jected data using the symmetries of the axisymmetric 
system so that all stars have positive coordinates. Then 
the data are again binned along the projected radius R. 
The measured velocity dispersion profiles are plotted as 
dots as a function of the projected radius R in Figure [H 
The sampling errors are not shown because they are very 
small. 

The velocity moments are related to the underlying 
potential <& and the distribution of the tracer v via a 
set of Jeans equations which in general read (Binney & 
Tremaine 2008) 

7nr(»4) + v^r + ™i*z- = -"inr- ( 4 ) 



dxi 



at 



dxi 



dxj 



Together with the continuity equation these equations 
form a set of four equations relating nine unknown func- 
tions, the three components of the mean velocity vj and 
six components of the symmetric tensor of- . The system 
can only be closed if some additional assumptions are 
adopted, as in the well-studied cases of spherical symme- 
try and axisymmetry, which we discuss below. Another 
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Fig. 10. — The velocity dispersion profiles of the stellar compo- 
nent measured in spherical coordinates r, 9, <j> in shells containing 
equal number of stars for models A (upper left panel), B (lower left 
panel) and C (lower right panel). In each of these panels the solid 
line refers to a r , the dashed line to ag and the dotted line to a^. 
The upper right panel shows the velocity dispersion profiles of the 
stellar component in dwarf A in cylindrical coordinates R, <fi and z 
as a function of R. The solid, dashed and dotted line refers to an, 
a^ and a z respectively. 
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Fig. 11. — Left panel: the spherical anisotropy profiles /3 of the 
three simulated dwarfs A, B and C (thin solid, dashed and dotted 
line respectively) and the axisymmetric anisotropy parameter k 
(reduced by unity, thick solid line) for dwarf A. Right panel: the 
relative differences between dispersion profiles (a^ — ag)/ag for 
dwarfs A, B and C (thin solid, dashed and dotted line respectively) 
and (an — a z )/a z (thick solid line) for dwarf A. 



possibility, for which Jeans equations are solvable, is the 
case of Stackel potentials (van de Ven et al. 2003). We 
find however that none of the well-known triaxial po- 
tentials of this form (de Zeeuw et al. 1986) reproduces 
the density distribution of all our dwarfs. In particular, 
none of the known models reproduces the rather com- 
plicated structure of dwarf A, although the distribution 
of mass in dwarfs B and C is well fitted by the perfect 
prolate spheroid (a variation of the so-called perfect el- 
lipsoid, see de Zeeuw 1985). We will further explore the 
similarities of those dwarfs to the perfect ellipsoid mod- 
els elsewhere (Lokas et al., in preparation). In any case, 
with the kind of kinematic samples available for dSph 
galaxies at present or in the near future there is little 



hope of constraining any models beyond spherical or ax- 
isymmetric (Statler 1994a, b; van den Bosch & van de 
Ven 2009). 

6.2. Spherical models 

We first model only the velocity dispersion profiles us- 
ing spherical models. For this purpose we solve the lowest 
order Jeans equation for spherical systems 



— {pa r ) + —pa r 



d$ 

^ dr 



(5) 



where we will assume that mass follows light so v = p 
up to a constant factor. The radial velocity dispersion 
ay is related to the angular dispersions ag and a$ by the 
anisotropy parameter /3 = 1 — u 2 , jo 2 . For simplicity we 
will consider only models with (3 =const i.e. given by the 
distribution function f(E,L) = f{E)L~ 213 . Such models 
require that there is no rotation and ag = a^. The ac- 
tual radial dependence of the dispersions (measured with 
respect to the mean) oy , ag and a^ on radius is shown in 
Figure [10] for our dwarfs A, B and C. The corresponding 
profiles of (3, calculated by replacing cr 2 with (ag+a 2 ) )/2, 
are plotted in the left panel of Figure [TT1 The right panel 
of the Figure shows the relative differences (er^ — ag)/ag. 
We can see that the condition ag = a^ is best obeyed 
for model B and although f3 profiles vary slowly with ra- 
dius, our simulated dwarfs are not strongly anisotropic. 
We believe this level of anisotropy should have negligible 
effect on our findings. 

For comparison with observations the solutions of 
equation ((5|) have to be projected along the line of sight 
(see e.g. Lokas 2002; Lokas et al. 2005). As demon- 
strated in Figure [21 our models have mass-to-light den- 
sity ratio almost constant with radius so in the model- 
ing we assume that the mass follows the distribution of 
stars. Note that if the dwarfs were not strongly affected 
by tides or were in isolation, more general models with 
extended dark matter profiles would be more appropri- 
ate (e.g. Lokas 2002; Koch et al. 2007a; Battaglia et al. 
2008; Walker et al. 2009c). The 3D distribution of stars 
can be obtained by deprojection of their surface density 
distribution. The surface density profiles measured for 
the observations along the x, y and z axis for models 
A, B and C are plotted as dots in Figure [i"2l We have 
also shown that the density distribution of stars is well 
approximated by the Plummcr law (JTJ). The projected 
density distribution then has a very simple form 



E(i2) 



p(r) r 



R 2 



:dr = 



a 2 M 



n{a 2 + R 2 f 



(6) 



The fits of this formula to the surface density distribu- 
tions measured for the three lines of sight are plotted in 
Figure [T2l as solid lines and the fitted scale lengths a arc 
given in each panel of the Figure and in the first nine 
rows of Table [5] (We actually adjusted 2 parameters, 
the scale length and the normalization, e.g. the stellar 
mass, but the latter is not relevant here although it was 
used in the determination of the half-light radii listed 
in Table [H The fits were done with the same weighing 
scheme as before, only now applied to surface densities.) 

Table [5] summarizes the results of fitting the solutions 
of the Jeans equation ([5]) (by the standard x 2 minimiza- 
tion) to the velocity dispersion profiles seen along the x, 
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TABLE 5 

Fitted parameters of the surface density profiles. 



Fig. 12. — Surface density profiles of the stellar component observed along the x, y and z axis (respectively from the left to the right 
column) for dwarfs formed in simulation A, B and C (from the upper to the lower row). The filled circles show the measurements made 
along the projected radius R. The solid lines show the best fits obtained with the Plummer surface density profile J6} with parameters 
listed in Table [5] for the spherical case. 

ulated dwarfs within this radius. The last column lists 
the values of the goodness-of-fit measure \ 2 /N rescaled 
to values corresponding to the total sample of a thousand 
stars and a hundred stars per bin. Next, keeping the a 
value fixed we adjust M and (3 as free parameters. The 
corresponding solutions are plotted in Figure [5] as solid 
lines. A third option is to keep (3 = and adjust M and 
the scale length a. The resulting solutions are plotted 
in Figure [9] as dashed lines. We can see that in most 
cases by including one more parameter beside the mass 
we can significantly improve the quality of the fits. Note 
that (3 and a cannot both be constrained because they 
are strongly degenerate with each other, but any of these 
parameters improves the fits since they both control the 
shape of the dispersion profiles while the mass acts as a 
normalization. 

The mass of dwarf A is underestimated in all ap- 
proaches, the least so when the line of sight is along 
the longest axis. This is obviously due to neglecting the 
rotation and fitting only the dispersion profiles. This 
problem can only be self-consistently solved by consid- 
ering axisymmetric models with intrinsic rotation, as we 
do in the next subsection. The best-fitting anisotropy is 
zero to mildly radial in rough agreement with the aver- 
age value measured in 3D for this dwarf (3 = 0.22 ± 0.16 
where the uncertainty is the dispersion of [3 values in 



simulation 


model 


line of sight 


a [kpc] 


d[kpc] 


A 


spherical 


X 


0.55 








v 


0.76 








z 


0.87 




B 


spherical 


X 


0.46 








y 


0.73 








z 


0.74 




C 


spherical 


X 


0.72 








V 


0.96 








z 


0.96 




A 


axisymmetric 


X 


0.56 


0.77 






y 


0.89 


1.89 



y and z axis for models A, B and C. In each row of the 
Table the values of the parameters that are kept fixed are 
given in boldface. We adopted three approaches. First, 
assuming the shape of the mass profile given by the fit- 
ted scale lengths a we adjusted the solutions of the Jeans 
equation for the isotropic case (with the anisotropy pa- 
rameter (3 = 0) i.e. fitting only the total mass M. The 
corresponding solutions are plotted as dotted lines in Fig- 
ure O The sixth column of Table [5] gives the masses 
within r < r max divided by the actual mass of our sim- 
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TABLE 6 

Fitted parameters of the spherical models. 



simulation line of sight 


a [kpc] 


P 


M[10 7 M Q ] 


Jft™. VmaxJ 


X 2 /N 


A x 


0.55 





3.80 


0.91 


0.90 




0.55 


0.0020 


3.80 


0.91 


1.02 




0.70 





4.41 


0.99 


0.84 


V 


0.76 





3.25 


0.71 


0.45 




0.76 


-0.051 


3.25 


0.71 


0.48 




0.96 





3.76 


0.74 


0.35 


z 


0.87 





3.19 


0.66 


0.41 




0.87 


0.20 


3.26 


0.67 


0.05 




0.61 





2.67 


0.63 


0.03 


B x 


0.46 





2.35 


1.03 


0.94 




0.46 


0.26 


2.45 


1.07 


0.34 




0.29 





1.93 


0.90 


0.15 


y 


0.73 





2.11 


0.80 


0.63 




0.73 


0.22 


2.16 


0.82 


0.13 




0.48 





1.73 


0.75 


0.07 


z 


0.74 





2.12 


0.80 


1.12 




0.74 


0.30 


2.20 


0.83 


0.20 




0.42 





1.63 


0.73 


0.09 


C x 


0.72 





1.02 


1.56 


1.27 




0.72 


-0.20 


1.04 


1.59 


1.10 




1.65 





1.97 


1.20 


0.40 


y 


0.96 





0.92 


1.13 


0.40 




0.96 


-0.078 


0.92 


1.13 


0.39 




1.43 





1.23 


0.94 


0.24 


z 


0.96 





0.48 


0.59 


0.85 




0.96 


0.24 


0.47 


0.58 


0.31 




0.59 





0.36 


0.61 


0.34 



Note. — The goodness-of-fit measures \ /N are rescaled to values 
stars per bin. The fixed parameters are given in boldface. 

the left panel of Figure [TTJ For dwarf B the mass is typi- 
cally overestimated when the dwarf is observed along the 
longest axis and underestimated if seen along any other 
axis. Note that only in this case the fitted anisotropy 
is mildly radial for all lines of sight and in good agree- 
ment with the 3D value of = 0.33 ± 0.14. This is 
probably due to the fact that for this case the condi- 
tion ag = <70 required by the spherical models is best 
fulfilled (see Figure [10] and the right panel of Figure [TT|) . 
For dwarf C the mass is strongly overestimated when the 
dwarf is observed along the longest axis and underesti- 
mated along the shortest axis. The inferred anisotropy 
varies from mildly tangential to mildly radial, while the 
3D value is p = —0.08 ± 0.17. In this case the estimates 
are least reliable as expected from the highly perturbed 
state of this dwarf. Note that the best-fitting values of 
a differ significantly from the ones obtained by fitting 
the surface density profiles of the stars (Table [5|) . This 
approach would therefore only be justified if the density 
profiles were unknown or highly uncertain. However, the 
trends in mass estimates are similar to the case when a 
are fixed and anisotropy is fitted instead. With just three 
cases analyzed, out of which only one (dwarf B) is really 
reliable (dwarf A has significant rotation and dwarf C is 
close to disruption) we cannot draw any firm conclusions 
concerning the reliability of our estimates of p. However, 
the trend of mass estimates increasing with line of sight 
closer to the major axis is clear in all cases. 



corresponding to a total sample of a thousand stars with a hundred 

6.3. Axisymmetric models 

So far we have modeled only the line-of-sight velocity 
dispersion profiles assuming that the galaxy is spherical. 
We will now try to reproduce both the rotation veloc- 
ity profile and the velocity dispersion for our dwarf A 
and those lines of sight where the rotation was detected 
(observation along the x and y axis). Both moments 
are shown as a function of the distance along the ma- 
jor axis of the galaxy image in Figure [T3l For the fit- 
ting the mean velocities were assigned standard errors of 
the mean a/ \fn. For the purpose of this study we will 
assume that the galaxy can be approximated as an ax- 
isymmetric system and again that the mass distribution 
follows that of stars. Given that the flattening in our 
simulated dwarfs is much more pronounced in the stellar 
component (see Figure [T] and fig. 13 in Klimentowski et 
al. 2009a) than in the dark halo (which remains almost 
spherical: the densities measured along the major and 
the minor axis never differ by more than a factor of two) 
it does not seem to be a realistic assumption. However, 
such one-component models may still be able to repro- 
duce data better than spherical models. 

For the axisymmetric systems it is convenient to work 
in cylindrical coordinates (R, </>, z) where R is now the 
radial coordinate measured from the center of the galaxy 
in the equatorial plane. We will assume that the only 
streaming motion in the galaxy is the rotation so Wr = 
vZ = 0, 7^ 0. If we further assume that the system is 
described by a distribution function that depends only 
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Fig. 13. — The rotation curves (upper row) and the line-of-sight 
velocity dispersion profiles (lower row) for the observation along 
the x and y axis (left and right column respectively) for simulated 
dwarf A. The dotted lines show the best fitting solutions of the 
single component isotropic (k = 1) models with scale parameters 
a and d adopted from the surface density distribution of stars and 
only mass M fitted. The dashed lines correspond to best fitting 
models with a, d and M fitted. The solid lines are for the fits with 
M and k as free parameters. The best-fitting parameters are listed 
in Table [7] 




Fig. 14. — Surface density profiles of the stellar component of the 
dwarf in simulation A seen along the x and y axis (the left and the 
right panel, respectively). In both panels the filled circles show the 
measurements made in a strip of width ±0.2 kpc along the major 
axis of the image and the open circles the similar measurements 
along the minor axis. The solid and dashed lines show the best fits 
obtained with the projected density associated with potential ((3j 
with parameters listed in Table [5] for the axisymmetric case. 



on the energy and the z-component of the angular mo- 
mentum vector, / = f(E,L z ) then the second velocity 

moments obey v R = v%, the mixed moments vanish and 



the Jeans equations take the form (Binney 
2008) 



dR 



H dR 

d * 

oz 



Tremaine 

(7) 
(8) 



The rotation may be introduced in an arbitrary way, 
for example we may have (Satoh 1980; Binney et al. 



1990) 



k 2 {v 2 



(9) 



where k is the anisotropy parameter and k = 1 corre- 
sponds to the velocity distribution which is isotropic ev- 
erywhere in the system, namely v 2 R — v 2 — v 2 ^ — 1T^ 2 . 
With this assumption and the known potential and den- 



sity, equation ([8]) can be solved for 



.,2 



and intro- 



duced into equation ([7]) to obtain the rotation v$. For 
comparison with observations all quantities have to be 
projected along the line of sight (see Satoh 1980). 

The upper right panel of Figure [10] shows that our sim- 
ulated dwarf is not strongly anisotropic. In this Figure 
we plot the profiles of the three velocity dispersions mea- 
sured in cylindrical shells along R containing equal num- 
bers of stellar particles. We can see that they are very 
similar. The only significant difference is visible at R < 1 
kpc where <tr is somewhat larger than the other two pro- 
files. This is due to the radial motion of the stars in 
the inner remnant bar. The relative difference between 
the dispersions of velocities along R and z is further il- 
lustrated in the right panel of Figure [11] as the thick 
solid line. Its positive value means that the underlying 



condition of our axisymmetric models, 



,,2 — 



;? is not 



strictly obeyed. The profile of the anisotropy parameter 
k (reduced by unity), calculated by replacing v R with 

(v 2 R + v 2 )/2 in equation ©, is plotted as the thick solid 
line in the left panel of Figure [11] 

The density distribution in these axisymmetric systems 
will be described by the modified Plummer model 
with b — c = 

Z) = ~{(R 2 + z 2 + a 2 ) 2 + d 2 z 2 }^ (10) 

where R 2 = x 2 + y 2 , so the model will be fully deter- 
mined by its total mass M and two scale lengths, a and 
d. The ratio of the scale lengths d/a controls the amount 
of flattening and rotation in the system. We will again 
start with isotropic [k = 1) models where the only free 
parameter is the mass, while the scale-lengths are deter- 
mined from the surface density distribution of the stars. 
To get them we measure the surface density but this time 
separately along the major and minor axis of the galaxy 
image. The results for the two lines of sight are illus- 
trated in Figure rj3] These profiles were fitted with the 
projections of the density following from (|10[) and the re- 
sulting scale-lengths a and d are listed in the two lower 
rows of Table [5] The quality of the fit is much better for 
the observation along x (left panel of Figure [14]) than for 
the observation along y. This is not surprising since in 
the latter case the bar is more pronounced. 

We proceed to solve the Jeans equations for an axisym- 
metric system and fit the solutions to the rotation curves 
and velocity dispersion profiles measured along the x and 
y axis and plotted in Figure 1131 Table [7J summarizes the 
results. As for the spherical models, the fixed values of 
the parameters are given in boldface and the X 2 /N val- 
ues are rescaled to correspond to realistic samples. The 
Table lists values obtained for different assumed inclina- 
tions of the galaxy where i is the angle between the rota- 
tion plane and the line of sight. The best-fitting profiles 
of the velocity moments for i = 0° are plotted in Fig- 
ure 1131 Again, we tried three approaches, starting with 
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TABLE 7 

Fitted parameters of the axisymmetric models for simulation A. 



line of sight 


i[dcg] 


a [kpc] 


CL [KpCJ 


h. 


ivl [-LU IVlQj 




A/ 2 IN 


X 





0. 


.56 


0.77 


1 


4.42 


1.03 


0.74 







0. 


,56 


0.77 


1.08 


4.52 


1.05 


0.63 










.54 


0.76 


1 


4.31 


1.01 


0.81 




10 


0, 


.56 


0.77 


1 


4.43 


1.03 


0.76 




10 


0. 


,56 


0.77 


1.09 


4.54 


1.06 


0.61 




10 





.54 


0.77 


1 


4.31 


1.01 


0.81 




20 






0.77 


1 


4.46 


1.04 


0.82 




20 


0, 


,56 


0.77 


1.12 


4.60 


1.07 


0.56 




20 





.54 


0.79 


1 


4.35 


1.02 


0.82 




30 


0. 


,56 


0.77 


1 


4.50 


1.05 


1.00 




30 


0, 


.56 


0.77 


1.18 


4.70 


1.09 


0.49 




30 





.53 


0.84 


1 


4.44 


1.04 


0.84 




40 


0, 


.56 


0.77 


1 


4.57 


1.06 


1.45 




40 





,56 


0.77 


1.28 


4.85 


1.13 


0.40 




40 





.53 


0.90 


1 


4.60 


1.07 


0.90 




50 


0, 


.56 


0.77 


1 


4.62 


1.08 


2.39 




50 


0, 


.56 


0.77 


1.46 


5.04 


1.17 


0.32 




50 





.52 


0.98 


1 


4.85 


1.13 


1.07 


y 





0. 


.89 


1.89 


1 


4.38 


0.81 


3.85 







0, 


.89 


1.89 


0.84 


4.11 


0.76 


2.28 










.95 


1.39 


1 


4.58 


0.85 


0.62 




10 


0. 


,89 


1.89 


1 


4.42 


0.81 


3.76 




10 





,89 


1.89 


0.84 


4.14 


0.76 


2.25 




10 





.95 


1.40 


1 


4.60 


0.86 


0.63 




20 


0. 


,89 


1.89 


1 


4.52 


0.83 


3.44 




20 





,89 


1.89 


0.85 


4.25 


0.78 


2.15 




20 





.95 


1.45 


1 


4.65 


0.86 


0.64 




30 


0, 


.89 


1.89 


1 


4.70 


0.87 


2.84 




30 


0, 


.89 


1.89 


0.87 


4.46 


0.82 


1.95 




30 





.95 


1.53 


1 


4.78 


0.88 


0.68 




40 


0, 


.89 


1.89 


1 


4.98 


0.92 


1.98 




40 


0, 


.89 


1.89 


0.91 


4.81 


0.89 


1.65 




40 





.95 


1.67 


1 


5.02 


0.91 


0.76 




50 


0, 


.89 


1.89 


1 


5.36 


0.99 


1.20 




50 





,89 


1.89 


0.99 


5.34 


0.98 


1.26 




50 





.95 


1.85 


1 


5.45 


0.98 


0.96 



Note. — The goodncss-of-fit measures \ 2 /N are rescaled to values corresponding to a total sample of a thousand stars with a hundred 
stars per bin. The fixed parameters are given in boldface. 



fitting the mass only with a and d adopted from the fit 
to the surface density profile and k = 1. The best- fitting 
solutions are plotted in Figure [13] as dotted lines. 

The fits are quite good in the case of the observation 
along the x axis but rather poor in the case of observa- 
tion along the y axis, especially for the rotation curve. 
The rotation is much too large compared to the data in 
spite of fitting the mass. This overshooting of rotation 
is due to a very large d/a ratio coming from the fit of 
the surface density distribution of the stars. In other 
words, self-consistent axisymmetric and isotropic models 
would predict more rotation for these d/a ratios. Still, 
the masses within r max found in this case are reasonable 
(see Table [7]), i.e. they are underestimated by no more 
than 20 per cent. 

As for spherical models we next relax the assumption of 
isotropy and fit the data again with two free parameters, 
M and k. As discussed in detail by Binney et al. (1990), 
decreasing k can suppress the rotation curve of the mod- 
els while increasing the dispersion profile. The results of 
the fitting are plotted in Figure [TBI as solid lines. As ex- 
pected, the quality of the fits is now improved, especially 



for the observation along the y axis (right hand panels 
of Figure I13[) where the rotation curve is now in better 
agreement with the data. Note however, that the veloc- 
ity dispersion profile is now worse reproduced and also 
the masses within r max are even more underestimated in 
comparison with the true values. A better approach in 
this case may therefore be to keep k = 1 and adjust a 
and d instead together with the mass. The corresponding 
profiles of the velocity moments are plotted in Figure [TS] 
as dashed lines. The fits are now considerably improved 
and the masses are not underestimated by more than 15 
per cent. Note also that the fitted values of a and d are 
not very different from those estimated from fitting the 
surface density distribution of the stars. 

Figure [15] shows the surface density contours and the 
line-of-sight kinematics of these best-fitting axisymmet- 
ric models in a way analogous to Figures [5] and [8] We can 
see that the models qualitatively reproduce the shapes 
and kinematics of the simulated dwarf A. In particular, 
they reproduce the faster variation of the velocity mo- 
ments with projected distance when the dwarf is seen 
along the longest axis, compared to the perpendicular di- 
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Fig. 15. — The structure and kinematics of the best fitting ax- 
isymmetric models for the observation along the x (left column) 
and y axis (right column) with fitted parameters a, d and M listed 
in Table [7] assuming inclination i = 0° . The upper row panels 
show the contours of the surface density £ in Mq kpc -2 plotted 
at AlogS = 0.2 with the innermost contours corresponding to 
logE = 7.6 (left column) and 7.2 (right column). The middle and 
lower row panels show the contours of the line-of-sight velocity and 
velocity dispersion respectively, with contour levels in km s — 1 in- 
dicated in each panel. Negative velocity contours are plotted with 
dashed lines. 

rection. Obviously, they are unable to reproduce the ve- 
locity gradient seen along the projected minor axis in the 
middle panel of the upper row in Figure[5]since this would 
only be possible with triaxial models (Binney 1985). 



7. SUMMARY AND DISCUSSION 

Using high-resolution A^-body simulations we have 
studied the stellar properties of dSph galaxies resulting 
from the tidally induced morphological transformation 
of disky dwarfs on a cosmologically motivated eccentric 
orbit around the Milky Way. We find that although the 
simulations differed only by the initial inclinations of the 
disk with respect to the orbital plane (0°, 45°, and 90°), 
this single parameter change led to significantly dissim- 
ilar outcomes. The reason for this is the well-known ef- 
fect that the tidal force is more effective in removing 
stars which are on more prolate orbits with respect to 
the orbital motion (e.g., Read et al. 2006). The dwarf 
spheroidal galaxies obtained had different masses, den- 
sity distributions, and velocity structure. Although we 
have used a most probable orbit of dwarf galaxies in the 
Local Group (Diemand et al. 2007), the variation of or- 



bital parameters would further differentiate the outcome 
of the simulations. On more circular orbits the dwarfs 
would be less affected by impulsive heating at percen- 
ters and therefore we expect them to be less evolved, 
preserving a more prolate shape and more rotation as in 
earlier stages of the evolution on our eccentric orbit (see 
Klimentowski et al. 2009a). 

The least evolved dwarf A (with initial inclination of 
90°) has an intrinsically triaxial distribution in the stellar 
component and preserved a fair amount of rotation. The 
rotation is mostly around the shortest axis. The inter- 
mediate case, dwarf B (with initial inclination of 45°), 
is a prolate spheroid and retained little rotation. The 
shape of its stellar component is still rather regular and 
so the dwarf can be considered a relaxed object. The 
most strongly evolved dwarf C (with initial inclination 
of 0°) is also a prolate spheroid with very little rotation 
but its noisy density contours and profiles of velocity mo- 
ments indicative of small numbers of bound stars suggest 
that it is very close to disruption. 

Despite the aforementioned differences and the fact 
that the final masses of the simulated dwarfs differ signif- 
icantly (see Table Q] and Figure [5]), the final M/L ratios 
are quite similar and fairly constant with radius. The 
former is related to the fact that the initial mass-to-light 
distributions of the dwarfs were identical and indicates 
that the dark matter and stars are lost in proportional 
amounts in all cases. The latter is a consequence of both 
our choice of initial conditions for the dwarf galaxy mod- 
els and the fact that most of the dark matter outside of 
the stellar component is stripped efficiently. Indeed, in 
the present study we adopt the Mo et al. (1998) scal- 
ings which connect the angular momentum/size of the 
disk with those of the parent dark matter halo. These 
general scalings together with our specific choice of pa- 
rameters result in an initial dwarf galaxy with a M/L 
ratio which is constant within 2 kpc or so and increases 
thereafter (Klimentowski et al. 2007). As the outer dark 
matter halo gets tidally stripped, the effect of increasing 
M/L ratio at larger distances, characteristic of the ini- 
tial conditions, disappears. However, it is not difficult to 
imagine scenarios where this situation can change as, for 
example, in cases where the original disky dwarfs are con- 
structed with much smaller scale lengths compared to the 
scale radius of the parent halo, r s . Such modeling would 
probably result in a non-constant M/L ratio within the 
luminous radius of the dwarf as the outer dark matter 
halo is tidally truncated. Until cosmological simulations 
reveal the structure of the dwarf progenitors themselves, 
the choice of initial structural parameters will remain 
somewhat arbitrary. 

We studied the density distribution of the stellar com- 
ponent both radially averaged and along the principal 
axes. In the former case it is well fitted by the Plum- 
mer law, in agreement with real dSph galaxies. For the 
latter we proposed a simple formula for the potential of 
the stellar component of the simulated dwarfs based on 
the generalization of the Plummer law. This formula re- 
produces quite well the triaxial shape of model A and 
prolate spheroidal shapes of models B and C (with a 
smaller number of parameters in the latter cases) both 
when measured along the principal axes and in projec- 
tion. The surface density distributions are quite similar 
to the observed ones, although a detailed comparison is 
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beyond the scope of the present paper. 

We also investigated the velocity distribution of the 
stars in the dwarfs along the principal axes and along 
the line of sight. The velocity structure seems to be the 
richest in the case of model A due to its triaxiality. In 
this case we detect rotation both around the short and 
the long axis which manifests itself in the misalignment 
of the zero velocity curve with respect to the projected 
minor axis. The velocity structure is much simpler in 
dwarfs B and C since almost no net streaming motion 
is seen there. In all models the dispersion of velocity is 
largest for velocities along the major axis which is due to 
the origin of the shapes of the dwarfs. All of them went 
through a prolonged bar-like phase where radial orbits 
along the bar are dominant and were transformed into 
prolate spheroids by systematic shortening of the bar. 

The projected velocity dispersions of our simulated 
galaxies decline slowly with radius, although some ex- 
hibit larger gradients. This behavior agrees with the ob- 
served profiles in real dSph galaxies (e.g. Lokas et al. 
2008; Lokas 2009) when the samples are cleaned of con- 
tamination using algorithms such as the one described 
in Klimentowski et al. (2007). Note that the majority 
of velocity dispersion profiles published in the literature 
(e.g. Muhoz et al. 2005; Koch et al. 2007b; Walker et 
al. 2007) are rather flat, which in our opinion is proba- 
bly due to contamination from the Milky Way or tidally 
stripped stars but may also be caused by multiple stellar 
populations (Battaglia et al. 2006) or an extended dark 
matter component (e.g. Walker et al. 2006). 

This velocity structure has important consequences for 
mass modeling. We fitted the line-of-sight velocity dis- 
persion profiles of the dwarfs observed along the princi- 
pal axes using different sets of free parameters. When 
the observation is done along the longest axis (the ve- 
locity dispersion is then larger and declines more steeply 
with radius) the mass is typically mildly overestimated. 
When it is done along the intermediate or short axis, it 
is underestimated. 

The effect is largest for dwarf C, which is close 
to disruption, but even in this case the mass is 
over/underestimated by only about 60 percent. This sug- 
gests that the very high mass-to-light ratios estimated for 
some dSph galaxies, such as Draco or Ursa Major, can- 
not be explained by departures from sphericity, but are 
rather the result of the formation process of their progen- 
itors, or are caused by other mechanisms that affected 
their baryonic mass fraction such as photoevaporation of 
the gas after reionization (e.g. Babul & Rees 1992; Bul- 
lock et al. 2000) or the combined effect of heating by the 
UV background and subsequent ram pressure stripping 
in the hot gaseous corona of the Milky Way (Mayer et 
al. 2006; 2007). 

For dwarf B, which seems to provide the best model 
for the real dSph galaxies (it still has a regular struc- 
ture and shows little rotation) we find that spherical 
modelling may overestimate the mass by up to 7 per- 
cent (when the line of sight is along the major axis) but 
underestimate it by up to 27 percent (when viewed along 
a short axis). Large kinematic samples of a few hundred 
up to a few thousand stars with measured velocities re- 
cently obtained by Walker et al. (2009b) for Carina, 
Fornax, Sculptor and Sextans allow for mass determi- 
nations with statistical errors as small as 5-15 percent 



(Lokas 2009) using simple models (assuming constant 
mass-to-light ratios and anisotropy) described here or 
similar. This means that the systematic errors due to 
non-sphericity of dSph galaxies are comparable or larger 
than the statistical errors in mass determination due to 
sampling errors of velocity moments for the best-studied 
dwarfs. 

For dwarf A the velocity moments (the rotation curve 
and the velocity dispersion profile) are much better re- 
produced with axisymmetric models constructed under 
the assumption that only the rotation around the short- 
est axis in present (i.e. neglecting the rotation around 
the longest axis). Although the presence of rotation has 
been reported in a number of dSph galaxies, like Ursa 
Minor (Hargreaves et al. 1994; Armandroff et al. 1995), 
Leo I (Sohn et al. 2007), Sculptor (Battaglia et al. 2008), 
Carina (Muhoz et al. 2006) and recently for more distant 
ones, Cetus (Lewis et al. 2007) and Tucana (Fraternali 
et al. 2009), it remains to be seen if such models are 
applicable to real dSph galaxies as they would require 
measurements of rotation curves at higher significance 
level than presently available. 

We have demonstrated that the kinematics of our sim- 
ulated dwarfs depends sensitively on the line of sight. 
This has intriguing implications for the missing satellites 
problem (Klypin et al. 1999; Moore et al. 1999), which 
is formulated by mapping the observed stellar velocity 
dispersions of dwarf galaxies to subhalo maximum circu- 
lar velocities, V max (see also Klimentowski et al. 2009a). 
These results imply that because our remnant dwarfs are 
non-spherical, U ma x depends both on the specific model 
used to fit the potential of the dwarf and the line of 
sight used to measure stellar velocity dispersions, a. In 
our models, the measurement of the central velocity dis- 
persion alone, when done along the longest axis can be 
higher by 50 percent compared to the perpendicular di- 
rection. We note that these findings are only relevant 
to dSphs that have a moderate mass-to-light ratio like 
our simulated dwarfs (e.g. Fornax or Leo I). Dwarfs em- 
bedded in a much more massive dark matter halo may 
in principle have a different relation between a and V ma , x 
since the stars could probe a region well inside the radius 
at which V max occurs. 

Lastly, the findings of the present study demonstrate 
that for a given orbit and initial density distribution the 
stellar structure and kinematics of tidally stripped disky 
dwarfs depend sensitively on the initial inclination of the 
disk with respect to the orbital plane. Similar conclu- 
sions were reached by Mastropietro et al. (2005) using 
A-body simulations to follow the tidal evolution and ha- 
rassment of disk galaxies inside a galaxy cluster envi- 
ronment. The models presented in this paper as well as 
those of Mastropietro et al. (2005) do not include gas. 

Mayer et al. (2006, 2007) have shown that gas present 
in the disky progenitors of dSphs can be either com- 
pletely stripped (by ram pressure and tides) or partially 
retained producing periodic bursts of star formation at 
pericentric passages. The outcome depends sensitively 
on when the dwarf fell into the primary potential. The 
discriminating factors are the strength of the cosmic ul- 
traviolet background and the local flux from the primary 
at the time of infall. In cases when gas is partially re- 
tained, it can affect the structural properties of the rem- 
nants, favoring a longer lived stellar bar and reducing 
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tidal mass loss by increasing its central density following 
a bar-driven gas inflow. This would tend to produce a 
remnant more similar to that of model A of this paper, 
i.e. triaxial and with more residual rotation. However, as 
the gas is consumed and the dwarf becomes devoid of it, 
further tidal heating may change the internal structure 
of the dwarf again. Fornax, Carina and Leo I, that had 
prolonged star formation, may have gone through these 
different phases. Overall, the inclusion of gas is likely to 
induce further scatter in the final structural properties of 



dwarf galaxies orbiting the Milky Way. We will extend 
our analysis to gasdynamical models in future work. 
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